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We derive the two-loop effective action for covariantly constant field strength of pure Yang-Mills 
theory in the presence of an infrared scale. The computation is done in the framework of the 
worldline formalism, based on a generalization procedure of constructing multiloop effective actions 
in terms of the bosonic worldline path integral. The two-loop /3-function is correctly reproduced. 
This is the first derivation in the worldline formulation, and serves as a nontrivial check on the 
consistency of the multiloop generalization procedure in the worldline formalism. 
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I. INTRODUCTION 

The physics of strong but slowly varying chromomag- 
netic and electric fields may provide some insight to the 
non-trivial vacuum structure of QCD. For example, the 
Euler-Heisenberg Lagrangian in QED exhibits truly non- 
perturbative effects; it allows the investigation of the non- 
linear regime of QED, and has also been studied beyond 
one loop, see e.g. [J i ||. In QCD non-linearities al- 
ready are present on the classical level due to its non- 
Abelian nature. It is yet unknown how the full effective 
action changes beyond one loop. The computation of 
multiloop terms in the effective action is cumbersome, in 
particular for a non-Abelian gauge group. These com- 
putations simplify if background field methods are em- 
ployed and the multiloop terms are evaluated for specif- 
ically chosen background configurations such as covari- 
antly constant or selfdual configurations. This gives ac- 
cess to beta functions and parts of the effective action 
beyond one loop, see e.g. 3 , [ M @| . Wordline methods 
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been shown to lead to a striking simplification of certain 
computations. A summation of Fcynman diagrams is 
already implemented without loop momentum integrals 
and Dirac traces [H E| ■ 

In the present work we provide for the first time a nu- 
merically accessible expression for the two loop effective 
action of Yang-Mills theory for covariantly constant fields 
in the presence of a physical infrared cut-off (at finite 
correlation length). This computation serves many pur- 
poses: firstly, it completes the construction of multiloop 
worldline methods for Yang-Mills theories as initiated in 
[l3l [T^ |. The two-loop beta function serves as a non- 
trivial consistency check. Secondly, it is the necessary 
input for an RG- improved non-perturbative computation 
of the effective action within a Wilsonian framework. Fi- 
nally one can study the stabilisation of the Savvidy vac- 
uum beyond one loop. 



II. WORDLINE REPRESENTATION BEYOND 
ONE LOOP 

We briefly recapitulate the analysis of [H The 
starting point is the generating functional of pure Yang- 



Mills theory in the presence of a background field config- 
uration A, 



m " / 



Da exp — S[a, A] 



where the gauge-fixed Yang- Mills action is given by 
S[a,A] = X - ftrF^(a + A) 2 



(1) 



(2) 



Tr \n(-D{A)D{a + A)) , 



with txt a t h — —5 ab /2 in the fundamental representation, 
D(A) = d — igA, and gF^ = i[D^,D v ]. Within a per- 
turbative expansion Z\A\ reads 



Z[A]=Y d Z n [A] 



(3) 



where Z n comprises the n-loop contribution to the gener- 
ating functional. The generating functional ([3]) is gauge 
invariant under the gauge transformation A — > A + Dto, 
and its logarithm is the gauge invariant Wilsonian effec- 
tive action of pure Yang-Mills. The one loop contribu- 
tion Z\[A] for a subset of field configurations, e.g. co- 
variantly constant fieldstrength F, has been computed 
in jl2l . for related results within standard methods see 
[2l| . So far, a full computation of the two loop contri- 
bution Z 2 [A] is lacking. After some algebra, |T]) can be 
turned into a more convenient representation for Z 2 [A], 
see (l3| . However, in all representations the formal ex- 
pression (p]) suffers both from UV and IR divergences. 
In the present work we regularise and renormalise these 
divergences separately: for the UV divergences we em- 
ploy dimensional rcgularisation for analytic purposes and 
construct a gauge-invariant proper-time cutoff interest- 
ing for numerical work. The divergences are then cured 
by appropriate counter terms. Additionally we intro- 
duce a physical IR cut-off. IR divergences are absent, 
if putting the theory into a box of size L. Effectively 
this can be implemented by introducing gauge invariant 
masses m ~ 1/L to the propagating degrees of freedom 
of the theory. The latter also has the advantage of im- 
plementing a physical mass gap on the level of the Green 
functions. This offers a path towards a self-consistent 



investigation of QCD in the confining regime in an effec- 
tive field theory approach that is quite close to the fun- 
damental theory. We emphasise that the approach leads 
to a fully gauge invariant effective action. However, the 
choice of the gauge-fixing parameter is directly related 



to the choice of different physical boundary conditions 
on the surface of the box. 

Within this framework the renormalised two loop con- 
tribution Z 2 [A] is provided by [f3| 



Z 2 \A] = exp(- / (-- - - - D 



A*3 



■(s/i,ya)(-| 



XeXp{ ^^ 1 fo^ }eXp{_ 5 TrlnA_1} exp{Trln(^ 2 +m 2 



a=/3=0 



+ c.t. , 



(4) 



with antisymmetric tensor field a^ v — —a vfJ ,, 

(A-Xi = (^-J) ab + 2f abc a^, (5) 

and 

= i-9^D 2 + 2F^f abc + rn 2 + (1 j)D^D v ]-^ .(6) 

We emphasise again that m 2 serves a twofold though 
related purpose. Firstly it accounts for a possible non- 
perturbative mass-gap, secondly it mimics the implemen- 
tation of a finite volume. The space-time integration over 
yi , j/2 can be used for regularising the generating func- 
tional ([3]) by means of the dimensional regularisation 




d n y , with D = 4 - 2e . 



In Feynman gauge £ = 1, the expression simplifies as 
the last term in ([6]) drops out, and the non-trivial tensor 



structure disappears. In Q we have not specified the 
counter terms indicated by c.t., that shall be discussed 
later. Within the representation introduced above the 
effective action derived from the generating functional 
flU reads 

r Fgl U on ~l~ Tghost j (7) 

with 

Tghost = ~ I -^A^Trln^ 2 + ™ 2 - £>/3/)l/3=o 
+c.t. , (8) 

the contribution of purely gluonic loops, r g i uon , has been 
computed in [14j |. and the result is summarised in Ap- 
pendix|B] We complete the analysis of [lj] by computing 
the ghost contribution r g h os t as well as the renormalisa- 
tion insertions. To that end we turn ([5]) into Euclidean 
wordline integrals, and arrive at [l3j 



x(Ti)=y' 2 /•2(T 2 )=i/; / -w(T 3 )=a:(T 1 ) 



Tghost = -~ / f dT 1( lT 2 dT 3 I [Vx] Tl I [Dx]t 2 I [Vw]t 3 

Jyi,V2 Jt Jx(0)=yi Jx(0)=y 2 Jw(0)=x(T 2 ) 

x(Pe /o T3 {Pe f T i N {x))cj Dl>(y' 2 )(\ e ) fg (Pef° 2N ^r D b ;(y[)(X a ) b c\ yi=yi +c.t.,(Q) 

I 

with the abbreviations Eq. ([9]) stands for the fully renormalised ghost contribu- 

tion to the two loop effective action. We have employed 

/[DxIt^Ix] = / Vx e~ f° dT (\ ±2 ( T )+ n i 2 ) p]^ an additional ultraviolet regularisation in the proper- 

J time integrations which entails a gauge invariant momen- 

tum cut-off. Such a cut-off scheme is amiable to numeri- 
an d cal computation, whereas the dimensional regularisation 



1 

N(x) ab = -iA c u x^(\ c ) ab . (10) 



facilitates analytic computations. In the following we 
shall conveniently project onto either regularisation by 
simply switching off either the regularisation parameter 



M®1{x(t)] — 2i(F° v — 5 11V -A c q x r '){\ c ) ab , shall conveniently project onto either regularisation by 



2 



t — > or 6 — ► 0. We have not specified the counter terms 
indicated by c.t., which in general depend both on the 
proper-time regularisation via r and on the dimensional 
rcgularisation via e. The computations of these counter 
terms will be discussed in the next section. 

For explicit computations we employ pseudo-Abelian 
su(2) with constant field strength, 

A^x) = A„(x)n a , with A lh = )-x v T Vil , (11) 

where n a is a constant unit vector in color space with 
n a n a — 1. The gauge fields (fTTj) satisfy the Fock- 
Schwinger gauge x^A^ = 0. With (TTTj) we rewrite the 
Lorentz matrices M and N as 



N(x) 



(12) 



with 71 — n a X a . The computation of r g h os t is straight- 
forward but tedious 221. It results in 



1 



ghost 



2(4^ J T 



dT\dT 2 dT 3 



x e~ m T 2g host [Ti, T 2 , T 3 ;T}+ c.t. , (13) 



with T = T\ + T 2 + T3. The integrand Ighost of the 
proper-time integral is given by 



Rghost 



— tr (TcotT^G^ 1 cos 2J-T 2 ) 



V ~ T 3 



(14) 



tr (G" 1 T 2 cot TT\ cot TT 2 ) 

J-tr iTG" 1 sin2JT 2 ) + tr (j^G" 1 )) det2#, 
J3 ' 



where 



G = TcotTT 2 + J 7 cot TT X 



(15) 



R = 



T 2 



sin FTi sin TT 2 ( 1 + .FT3 (cot TT X + cot .FT2 ) ) 



The expression (IT31 with (fTi)) is numerically accessible, 
after the counter terms in (fT3"|) are specified. 



III. RENORMALISATION 

Now we discuss the UV subtraction terms hidden in 
the counter terms that render Z 2 , T 2 finite, and in par- 
ticular, (|13[) finite. Apart from applying a standard di- 
mensional regularisation convenient for analytic consid- 
erations, we have introduced a gauge invariant UV regu- 
larisation by cutting off the proper-time integrals in (fl~3"l) 
at a finite proper time r, > r. This translates via 
a Laplace transform into a gauge invariant momentum 
cut-off if the effective action is formulated in terms of 



momentum loops. Such a scheme makes numerical com- 
putations accessible where the dimensional regularisation 
only can be employed in exceptional cases. Indeed, a fully 
non-perturbative wordline formulation of quantum field 
theories would provide a tool for devicing gauge- invariant 
momentum cut-off schemes on the non-perturbative level 
which would be highly interesting. The ghost action © 
is then written as r g h os t = r g host,rcg + c.t. with 



1 



- ghost, rcg 



dTidT 2 dT 3 e 



-171 T 



xX ghost [Ti,T 2 ,r3;^]. (16) 



The regularised expression r g host,rcg diverges if the reg- 
ularisation parameters r, e are removed. Here we first 
concentrate on the proper-time regularisation with e = 0. 
Then the regularised expression r g h st,reg in (fl6|) diverges 
with powers of 1/t, more precisely with l/r n (lnr) m . 
Moreover, since we are dealing with the two-loop effec- 
tive action, the divergent terms are not necessarily poly- 
nomial, and the counter terms cannot be determined in 
a polynomial expansion. The non-polynomial terms can 
be attributed to the divergence of one loop sub-diagrams 
which can be used to construct the related counter terms. 
However, here we want to set-up a procedure with which 
these counter terms can be derived systematically from 
r g host,reg by means of derivatives. Such a procedure mim- 
ics the standard BPHZ-type subtraction schemes in mo- 
mentum space. The divergences in r are extracted by 
appropriate r-derivatives with the help of the identity 

T^-T~ n {\nT) m = -r-"(lnr) m (n - . (17) 

Applying the above r-derivative to r g h os t,regi the physi- 
cally finite term drops out. Hence appropriate subtrac- 
tions 



r = r I cg-^c„(r)(r9 T )"r r c 



(18) 



render the action T finite. Note that such a procedure 
in principle only properly provides the renormalised ef- 
fective action by a careful discussion of the finite renor- 
malisation that originates in the subtractions in (fT8"|) . In 
the present two loop case it suffices to only take one r- 
derivative, e.g. 



Td T T, 



host, rcg 



1 dT.dT.dT^m-r) 
xe-^JghctlTi, T 2) T 3 ;^. (19) 



This reduces the number of proper-time integrations and 
makes the divergence structure analytically accessible. 
This procedure deserves further studies. 

We still have to compute the traces over the field 
strength. Following (25j , we work in the Lorentz frame in 
which the electric and magnetic fields are parallel, and 
thereby the field strength takes on a simple form with 



3 



only two non-zero symplectic block elements, i.e., it can 
be written as 



T = ao\ + ba 2 



with 



and 



a 





ai = 1 J ' CT2 = I a 



-IT], 



(20) 



(21) 



(22) 



where e and rj are the magnitudes of the magnetic and 
electric fields, respectively. 

We close with a remark on the explicit computation 
of the proper-time integrals. In particular for numeri- 
cal purposes it is advantageous to convert the integrals 
into less divergent expressions. Indeed, I g host as well 
as the corresponding integrand X g i UO n can be integrated 
analytically over T3 and hence can be written as a to- 
tal T^-derivative. The computations are deferred to Ap- 
pendix [X] and Appendix [B] respectively and the results 
read 



1 



j.Ik >sl 



2(4vr) 



D 



(23) 



x&^ighostPi, T 2 , T 3 ;T}+ c.t. . 
where Xghost is given in l|A.4[) , and 



2(4?r) 



1 r 00 _ 

-°l W Tie ~ m2T 



(24) 



xd T3 1 s Won[Ti, T 2 , T 3 ; JT] + c.t. 



where 2gi UO n is given in (|B. 1 1|) . The counter terms in 
([23| , (f24|) are r-dependent and can be computed from 
the derivative procedure outlined in (fT7|) . (TTS]). 



IV. TWO LOOP P FUNCTION 

In the remainder of this work we concentrate on the 
the question of full two-loop consistency of the wordlinc 
formalism suggested in [l3|, [3] , the construction of which 
we have completed here. To that end we discuss the run- 
ning of the coupling at two loop which is universal in 
mass-independent renormalisation schemes. As this con- 
cerns an analytic computation we employ a dimensional 
regularisation with r = 0. Moreover, we use a minimal 
subtraction scheme that renders the renormalisation con- 
stants mass- independent, and hence projects onto the the 
universal result for the two loop /3-function. How such 
a mass-independent scheme is fixed in the presence of 
general IR cutoffs has been discussed in detail in [23|, 
and we can straightforwardly use the related arguments 
for the mass terms for gluon and ghosts employed in the 



present work. The /3-function can be read-off from the 
running of the wave function renormalisation Z4 of the 
gauge field. Expanding the two loop contribution to the 
effective action I^rog in powers of T we are led to 



y 2 \a\ = ^- [ tr^ 2 + 0(^ 3 ). 



(25) 



The .^-coefficient of the ghost effective action r g host.rc 
reads 

4.9 4 , 10 , or , 8 3 



(47T) 



where the coefficients C[ are given by [l 
C( = (4V) 2e / ]JdT, 



' T^T 2 , T 



i=l 
3 



i=l 
3 



C 3 = (4V) 2e / WdT, 
Jo i=i 

poo 3 
J o i=l 



n 4 - 

T 1 T 2 T 3 -m 2 T 



, (27) 



with O = T X T 2 + T 2 T 3 + T 3 T X . The divergent parts of the 
coefficients C' read 



c 3 



0(e°), 



Ci = 



1 


v 9 3 




1 


/ 1 PrnA 

l 18 3 j e 


6? + 


1 


v 72 6 ; 


12e 2 


1 

127 + 


-O(e ), 



with 



Pm = IE + hi 



47T^ 2 



(28) 



7_e = Euler const . (29) 



Now we are in the position to perform the UV renormal- 
isation for the ghost term. Note that by power counting, 
the UV divergence can appear at most in the quadratic 
term in the expansion, we write 



Tghost,ren — r g i los t ^Cghost / ' 1 ^ 
+ ^Cfinitc / tr-^ 2 j 



(30) 



where 



_ ^9 f 10 , , 8 , 3 , 

°ghost - JfaMX 3"°1 ~ /L/ 2 - 3^3 + 2 4 



i\c[+ l -C' 2 + b -C', + C' i )e} 



(31) 
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and Cg nite is its finite part. This introduces the minimal 
subtraction scheme in the ghost part. The integral in 
(f!3| is changed by the additional integrands proportional 
to C' ghoBt J- 2 rendering a finite expression. We remark 
that it can be explicitly checked that the renormalisation 
constants are mass-independent at two loop. This con- 
stitutes a mass-independent RG-scheme and hence 02 is 
universal. 

The two loop /3-function is provided by 



= -g [id^lnZg 



(32) 



-9 ( PiCa [f n 



Wl(£)\o {9 «) 



The background field formalism allows us to directly ex- 
tract the two loop /3-function from Za'- the effective ac- 
tion T[A] is gauge invariant and consequently the combi- 
nation gA is RG-invariant, leading to Z g = 
hence 



-1/2 



and 



= A 1 df, In Z A . 
With ([SID and J33D we conclude that 



Z A = 1 



02 

e 



cl 



(33) 



+ 0(g 6 ), (34) 



and we directly read off the two loop /3-function from the 
subtraction terms computed in the last section. Using 
(|28p in (j2"6")l we arrive at the ghost contribution 02, ghost 
to the two loop coefficient 02, 







2, ghost 



ft 
IT 



(35) 



The gluon loop contribution has been computed in [14j 
as 02, gluon = —11/2. It is left to compute the contribu- 
tions of the one loop counter terms. They arise from the 
insertion of the one loop RG constants of coupling and 
propagating field at one loop. Note that the propagating 
field is the fluctuation field a with one loop wave function 
renormalisation 



SZ a = 



1 9 2 



e (4tt) 2 



1 



+ o(l"0 



(36) 



The worldline counter terms reduce to the standard one 
loop graphs. The corresponding diagrams shown in Fig.l 
and 2 result from the one-loop renormalisation of the 
fluctuation field a and its vertices, while those in Fig. 3 
and 4 arise from the one-loop renormalisation of ghost 
field and its vertices. 

The computation of the counter term in Fig.l requires 



a gluon mass renormalisation with m 
Feynman gauge it is given by [24| 



Z a Z m 



1 



i g 2 C A 

e (An) 2 



Z a Z m m. In 



(37) 



This counter term has been considered in [14j , which 
gives a contribution of 10/3e to the two loop coefficient 




FIG. 1: counter terms from the propagator renormalisation 
of the fluctuation field a . 




FIG. 2: counter terms from the vertex renormalisation . 



02- We also remark that the gauge- fixing term does not 
renormalise, £ — > Z a ^. 

Note that the ghost mass term is not renormalised. The 
counter terms from Fig.l-Fig.4 give rise to pole contri- 
butions 



Fig.l + Fig.2 : 
Fig.3 + Fig.4 : 



35 
1 

6i' 



(38) 



02 



(39) 



and we are led to 

11 11 35 10 1 

~T + lT + ir + T + 6 

which agrees with the well-known result, e.g. US- As we 
are only interested in diagrams with external background 
fields we could have rescaled the fluctuation field a and 



17 

y 



-□ 



V/VW^ 



FIG. 3: counter terms from the ghost propagator renormali- 
sation . 



FIG. 4: counter terms from the vertex renormalisation . 
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the ghost field with the renormalisation factors. With 
these rescaled fields, the diagrams above reduce to terms 
proportional to the renormalisation of the gauge-fixing 
term introduced by this rescaling and the renormalisa- 
tion of the mass terms which also changes by this rescal- 
ing (e.g. the ghost mass renormalises with these rescaled 
fields). This has been used in ||. Of course, this does not 
change the result. For comparison we list the different 
contributions 



/^-contributions 




m 2 ^0 




m 2 = 


two-loop diagrams 




n 

3c 




7 

3c 


Fig.1+2 


5 
18c 


, 80 35 , 
T 9e u 6e ~ 


10 
3 c 


10 
3 c 


Fig.3+4 


+ 


<r or 7^ 


+ 





total 




17 

3c 




17 

3c 



In the middle column the different contributions from 
a direct computation (right), and from one with rescaled 
fluctuation fields (left) are listed. 

We can use the above results on the consistent renor- 
malisation in the presence of an infrared mass-scale to 
define the renormalised two loop contribution r 2 [A] by 
means of appropriate subtractions instead of the di- 
mensional regularisation used above. This allows us to 
numerically compute the full two loop effective action 
r = Ti + r 2 as a function of T ' . 



V. OUTLOOK 

In the present work we have completed the wordlinc 
construction of the two loop effective action initiated in 
EH • In particular we have provided a crucial consis- 
tency check of the construction by computing the univer- 
sal two loop /3-function within the wordlinc formalism. 

We also have set-up a practical ultraviolet BPHZ-type 
renormalisation scheme in the proper-time which makes 
numerical computations accessible. For example, this can 
be used to numerically compute the two loop effective 
action for covariantly constant fields. 

The inclusion of fermions in the present approach is 
straightforward, and is, in our opinion, the physically 
most interesting extension of the present work. 
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A. GHOST CONTRIBUTION 

One of the T-integrations in (|13[) can be done ana- 
lytically. This is achieved by integrating T g host over T 3 . 
Performing all the traces in the integrand of (fT4]) and 
summing over them gives 



A 



2a 



-aT 3 + 2 cosh(2aT 2 ) coth(aTi) 



1 + aa\T 3 

+ aT 3 coth(aTi) coth(aT 2 ) + 2 sinh(2aT 2 ) 
+ {a <->&), (A.l) 

with a\ = coth(aTi) + coth(aT 2 ). The square-rooted 
determinant term reads 



B 



C'i 



(l + aa 1 T 3 )(l + 66iT 3 ) 



(A.2) 



where C x = a 2 b 2 csch(aTi) csch(&Xi) csch(aT 2 ) csch(&T 2 ). 
We define 

/•T3 

Tghost [Ti , T 2 , T3 ; T] — f dT 3 T g host [Ti , T 2 , T 3 ; T\ , 
Jo 

(A.3) 

after analytically performing the integration over T 3 one 
finds (bi = a\{a <-» b)) 



tghost 



Ad (ln(l + aai T 3 ) - ln(l + M»iT 3 )) 



ddi — bbi 



CiC 2 



>aai(aai — 6(?i) 2 (l + aa\T 3 ) 
- aa\{l + bbiT 3 ) ln(l + aaiT 3 ) 
+ aa 1 (l + 66iT 3 )ln(l + »iT 3 )) 

+ (a «+&))] -[T 3 ^0]}, 



( — aai + bbi 



(A.4) 



with 

C 2 = 2a 2 ( - 1 - 2ai cosh(2aT 2 ) coth(aTi) 

+ coth(aTi)coth(aT 2 ) - 2aisinh(2aT 2 )) . (A.5) 

B. GLUONIC CONTRIBUTION 

The gluon loop contribution to the two loop effective 
action reads I14i 



Tgluon = -^(4tt) 



dT x dT^dT 3 e~ m T T gluon + c.t 



(4n)- D I dT x dT^dT 3 e- mAT I det*(^— ) 



tr 



2 sin TT X cos 2T(T X + 2T 2 



6 



-2sin F(T X + T 2 ) cos^(2Ti + 3T 2 



-2 cos 2 T(T X + T 2 ) sin{TT 2 ) cos T(T± - T 2 ) 
T 

+ — [4sin^"Ti sin^Ta sin2.F(Ti + T 2 ) 

-2smTT x cos^(2Ti + 3T 2 ) 
-2sin^T 2 cosJ^(Ti - 2T 2 ) 

- sin^(Ti + T 2 ) cos 2^(Ti - T 2 ) 



r ^ 2 T, / 

+tr cos 2^T 2 tr ^— sin T{T X + T 2 ) 

lAjrsmj 7 l 2 \ 

x cos.F(2Ti + T 2 ) - smFTx cos2.F(Ti + T 2 )) 

+- r -(zwa.J : Tx cosJf(2Ti +T 2 ) 

+ cos2^Ti sin^(Ti + T 2 ) 
T 2 T, 



-tr 



-Ajr sin :FT 2 
-l-cos^Ta sin^XTi + T 2 ) 



sin-T-T^ cos .F(Ti — T 2 ) 



J" 



■ sin .FTi cos .FT 2 ) 



-sinJ^Ti cos2.fT 2 

xtr(cos2Jf(Ti + T 2 )j + 5(T 2 )2(1 - D)tr(cos 2TT X 
+J(T 3 )tr(cos2.F(Ti-T 2 )) 
-<5(T 3 )tr(cos 2JTi) tr (cos 2JT 2 

dTidT 2 / dets 



1 



+ 2^) 



— D 



JF 2 



sinJFTi sin TT 2 
ftrcos2J 7 (T 1 +T 2 ) - trcos2J 7 (T 1 - T 2 )) + c.t. (B.6) 



with 



A^ = an(JTi) sin(^T 2 ) + TT 3 sin[JT(T 1 + T 2 )] . (B.7) 



Following the same procedure as we extracted Eq. (I26| 
from (|13p. we obtain the renormalisation part of the ef- 
fective action above at the second order of T 



(i$ ( -i ) / ( -i v - )+o(e0) ' (R8) 



and hence the contribution of the gluon loops to the two 
loop /3-function coefficient is — ll/2e. 

For a numerical evaluation of r g i uon a less singular rep- 
resentation is advantageous. To that end we write r g i uon 



as 



gluon 



2(4*)" Jr i\ 



(B.9) 



'xd T3 I gluon [T 1 , T 2 , T 3 ; J 7 ] + c.t. 



with <9T 3 igiuon = 2giuou, and hence similarly to the ghost- 
contribution it reads 



Rgiuon Pi j T 2 j T3 j <F\ — I dT 3 X g i uon [Ti , T 2 , T 3 ; J 7 ] , 
Jo 

(B.10) 

The T3-integral in the definition of Zgi uon can be per- 
formed analytically and yields 



Rgiuon 



A'd (ln(l + aai T 3 ) - ln(l + M>iT 3 )) 



(CiC 3 0(T 3 ) 



2aCi(aiC 5 - C 4 ) 



66i V 1 J v J ' ai(aai — 6fei) 2 (f + 001X3) 

x ( - aai + 66i - aai(l + 6&iT 3 ) ln(l + aaiT 3 ) + ooi(l + 66iT 3 ) ln(l + &&1T3)) + (a <-> 6)) 



-[T 3 



0]}, 
(B.ll) 



with #(T 3 ) the step function and 



A' = 2a ( C 5 + C 4 aT 3 ) _ + CQsh ( 2ari ) 6 , T s + 2 g / r n / cos h(2a(Ti-r 2 ))-2(cosh(2aT 1 ) 
1 + aaiT 3 V 



+ cosh(26Ti)) cosh(2aT 2 ) - 2 sinh(2aTi) sinh(2aT 2 ) ) + ( a 4-+ b) . (B.12) 

The abbreviations C 3 , C4 and C5 read 

C 3 = 2 cosh (2a(Ti - T 2 )) - 4 sinh(2aTi) sinh(2aT 2 ) - 4( cosh(2aTi) + cosh(2&Ti)) cosh(2aT 2 ) , 
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C 4 = 4 cosh (26(Ti + T 2 )) ( - 1 + coth(oTx) coth(aT 2 )) + (2 cosh(26T 2 ) cosh (a(T x + T 2 )) 
+ 2 cosh (a(Ti - T 2 )) + cosh (o(3Tx + T 2 ))) csch(aTi)csch(aT 2 ) , 

C* 5 = 2 cosh (2a(Ti + T 2 )) coth(aT 2 ) + 2 cosh (26(Ti + T 2 )) coth(aT 2 ) - 2 cosh (a(Ti - 2T 2 ))csch(aTi) 
+ 6cosh(a(2r! + T 2 ))csch(aT 2 ) ( cosh(2aT 2 ) +cosh(2foT 2 )) - 2 cosh (a(2Tx + 3T 2 ))csch(aT 2 ) 

+ csch(aTi)csch(aT 2 ) sinh (a(Ti + T 2 ) (2 cosh(2aTi) cosh(2aT 2 ) + 2 cosh(2aTi) cosh(26T 2 ) 
- cosh (2a(Ti - T 2 ))) - 4sinh (2a(T x + T 2 )) . (B.13) 
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